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We investigate spherically symmetric self-similar solutions in Brans-Dicke theory. Assuming a 
perfect fluid with the equation of state p = (7 — (1 < 7 < 2), we show that there are no non- 
trivial solutions which approach asymptotically to the flat Friedmann-Robertson-Walker spacetime 
if the energy density is positive. This result suggests that primordial black holes in Brans-Dicke 
theory cannot grow at the same rate as the size of the cosmological particle horizon. 
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I. INTRODUCTION 

It has been pointed out that black holes could be 
formed in the early universe as a result of the collapse 
of the initial density perturbations, bubble collisions, or 
other mechanisms [Qj^,^- They are called primordial 
black holes (PBHs) and are of special interest since they 
could be the promissing candidates of MACHO or the 
only ones which could be small enough to be evaporating 
by the present epoch due to the quantum effects Q . The 
evolution of PBHs has been studied by many authors in 
general relativity 

Zeldovich and Novikov pointed out the catastrophic 
growth of the PBHs in a simple Newtonian argument 
which implies that the PBHs of which size is compara- 
ble to the cosmological particle horizon could continue 
to grow at the same rate as it throughout the radiation 
era, while those of which size is much smaller than that 
could not grow much at all [Q . One may think that the 
evolution of such PBHs could be described by self-similar 
solutions in general relativity. Carr and Hawking, how- 
ever, showed that there is no non-trivial spherically sym- 
metric self-similar solution approaching asymptotically 
to the flat Friedmann-Robertson- Walker (FRW) space- 
time at large radius for a radiation fluid This ar- 
gument was extended to the case of a general fluid with 
p = (7 - 1)^ (1 < 7 < 2) ||]. The result implies that the 
PBHs must soon become much smaller than the size of 
the cosmological particle horizon and cannot grow very 
much at all. 

It has also been argued that in the early universe, grav- 
ity might obey a scalar-tensor type theory rather than 
general relativity. A scalar-tensor type theory of gravity 
arises naturally as a low-energy limit of string theory. In 
such theories, gravitational "constant" is given by a func- 
tion of a scalar field which couples non-minimally with 
gravity so that the effective gravitational "constant" G 
generally varies in space or in time. Brans-Dicke the- 
ory is the simplest but still very important one among 



the scalar-tensor theories. Then, the evolution of PBHs 
in Brans-Dicke theories has been studied by several au- 
thors |,0,|ll|,[2|,[|,^]. 

Barrow proposed two scenarios for the evolution of the 
PBHs in Brans-Dicke theory Q . In his first scenario (A) , 
the Brans-Dicke scalar field evolves everywhere homoge- 
neously so that the gravitational "constant" at the black 
hole event horizon always the same as that at the cosmo- 
logical particle horizon. No information about formation 
epoch remains. In his second scenario (B), however, the 
gravitational "constant" at the black hole event horizon 
remains constant with time and is therefore always the 
same as that at the formation epoch, while it changes 
at distant scales. Since the gravitational "constant" at 
the horizon stays an old value when a black hole was 
formed, the remaining gravitational "constant" near a 
PBH is called the gravitational memory. Discussing the 
energy emission due to the Hawking radiation in these 
scenarios. Barrow and Carr found that there are the sig- 
nificant deviations from that in general relativity in both 
scenarios Two alternative scenarios, which are inter- 
mediate between the scenario (A) and (B), are proposed 
by Carr and Goymer | p^ . Analyzing perturbations of 
a scalar field in the Schwarzschild background jl^, Ja- 
cobson concluded that the gravitational memory is weak. 
However, this conclusion does not follow when the size 
of a black hole is comparable to that of the cosmolog- 
ical particle horizon. Using the Tolman-Bondi solution 
and neglecting back reaction of a scalar field to the back- 
ground geometry, Harada, Goymer and Carr studied the 
evolution of a scalar field when a black hole forms from 
the collapse of dust in the flat FRW spacetime and con- 



cluded that there is little gravitational memory |15|. 

The purpose of this short note is to study the evolution 
of the PBHs in Brans-Dicke theory by using self-similar 
solutions which contain a perfect fluid with the equation 
of state p = (7 — 1)^ (1 < 7 < 2). In this paper, we 
adopt the unit of c = 1. 
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II. SPHERICALLY SYMMETRIC SELF-SIMILAR 
SOLUTION IN BRANS-DICKE THEORY 

The field equations in Brans-Dicke theory are given by 
1 Stt 1 



2 \ ^'f^ 



T^"'.^ = 0, 



3 + 2lu 



(2.1) 
(2.2) 
(2.3) 



where the constant oj is the Brans-Dicke parameter and 
(j) is the Brans-Dicke scalar field, which is related to the 
gravitational "constant" by G = l/</). lo > —3/2 is re- 
quired for the theory not to contain tachyonic solutions. 

We consider a perfect fluid as a matter field T^i/ = 
PQtiv + + p)U^Uv, where is the 4- velocity. The ra- 
dial coordinate r is chosen to be comoving. The equation 
of state is assumed to be p = (7 — where 1 < 7 < 2. 
A solution of the Einstein equations is said to be self- 
similar if it admits a homothetic Killing vector ^ such 
that C^g^i, = 2g^i^, where denotes Lie derivative 
along ^. Cahill and Taub first investigated spher- 
ically symmetric self-similar solutions in which the ho- 
mothetic Killing vector is neither parallel nor orthogonal 
to the fluid flow vector. They showed that by the suitable 
coordinate transformations such solutions can be put in 
a form in which all dimensionless quantities are functions 
only of the dimensionless variable z = r/t. Generalizing 
their method to the present model, we find the follow- 
ing basic equations. The line element of a spherically 
symmetric self-similar spacetime is given by 

= -e2*(^)dt2 + e2*(-)dr2 + r''S''{z) dn\ (2.4) 

where dSl^ = dO^ -\- sir? Bdip^ . Here we assume the self- 
similarity of a perfect fluid and the Brans-Dicke scalar 
field independently, i.e. 



(2.5) 



where a, b and k are constants. These constants satisfies 
the relations a = b = —(2 — k ) th rough the equation of 
state and the field equations (pT|)-(2.5). Equation ( |2.5| ) 
implies that the quantities p, fi and (j> must be of the form 
as 



W{z) 



P 



P{z) 



(f>^r''tlj{z), 



(2.6) 



and then equations (2.1)-(p.3|) reduce to the ordinary dif- 
ferential equations w.r.t. the self-similar variable z. 

The energy-momentum conservation equation is writ- 
ten in self-similar variable as 



(2 - k)P - P' 7-1 



W + P 

' w 



W 



' s 



1 w 



(2.7) 
(2.8) 



where prime denotes the derivative with respect to Inz. 
These equations can be integrated to give 



(2-,t)(^-l) 



— cqz 1 'W 



(2.9) 
(2.10) 



wher e cp , ci are integration constants. Equations (2.1) 
and (2_^) then reduce to the following ordinary differen- 
tial equations; 

(2-hK)(7- 1) A„ , W'\ , 
- .) - — j - (. - 1 + u.) 



Ak — 1 + UK) — 



— — Ci5 W T, 



(2-^)(7-l) W'\ ( t. 

W J \ S V 



7 

„S" iw 

2 \ 

S -fW 



S' 

K + 2{ l + — 



0, 



(2.11) 



2^^ 

S V 



(2.12) 



-V + 2k - IH^ + K K - 1 



4 



{2-k){j-1) S' j~2W' 
S W 



1 



1 



(2-k)(7-1) -i-2W' 



ip ip \ 7 



7 w 



8^(37-4) ,^_4wi 



{3 + 2lo)^ 



(2.13) 



with V = ze*-*=(ci/co]£^-[(2-«)('r-i)/7] 5-2 pp^(7-2)/7. 
Equations (2.11) and (2.12) are independent of equations 
(P) and (pT 



HI. THE FLAT FRW SOLUTION 

The flat FRW solution in Brans-Dicke theory with a 
perfect fluid obeying the equation of state p = (7 — 1)/J, 
was found by Nariai JTzt and it is a particular solution 
of the reduced ordinary differential equations, i.e. a self- 
similar solution. The flat FRW solution in the self-similar 
variable {^FRwiz),^FRw{z), Sfrw{z),Wfrw{z),) is 



bo 



S 



FRW — 



Wfrw 



1 - q 

V;o(3 + 2a;)[2(5-37) + 3(2-7)^H ^3,7 
47rao2[4 - 37^(7 - 2)]2] ^ ' 



(3.1) 
(3.2) 

(3.3) 



2(4-37) 
4 - 37^(7 - 2) 



_ 2[l-u;(7-2)] 
'^-4-37^.(7-2)- ^^-^^ 
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One can put this solution in a more familiar form 

ds^ = -dt'^ + t'29(dr'2 + r'^dn^), (3.5) 

by making the coordinate transformation t' = aat, r' = 
a^'^bor^''^ /{I — q). In the case of radiation (7 — 4/3) , 



q — 1/2 and k = are obtained from equation (3.4). 
Since k = 2 — Sqj, k — implies that the Brans-Dicke 
scalar field is constant which case is equivalent to that 
in general relativity. Assuming the energy density is pos- 
itive (W > 0), we find 



LU > 



2(37 - 5) 
3(7-2)2 



(37 - 4)^ 
6(7 - 2)2 



(3.6) 



If 7 = 2, any value of lu is not possible. On the other 
hand, in the case of radiation (7 — 4/3), the energy den- 
sity is always positive for any value of cj > —3/2. 



IV. UNIQUENESS OF SELF-SIMILAR 
ASYMPTOTICALLY FLAT FRW SOLUTION 

We are interested in self-similar solutions which are 
asymptotically flat FR W so lution, i.e. in which W, S and 
ip approach the form (3J_)-(1_S) as z ^ 00. It is con- 
venient to introduce new functions A, B and C, which 
denote deviations from the flat FRW solution, defined 
by 



W ^WFBw{z)e^^'\ (4.1) 
S ^ SFRw{z)e''^'\ ^^^^FRw{z)e^^'^. (4.2) 

These functions denote deviations from the flat FRW 



solution, and equations (4.1) and (4.2) imply (through 
equations (2.7) and (^.8[)) that 



(4.3) 



The flat FRW solution is given by^ = S = C==0 



for any z. Since V = (60/00)2^ exp(— 2i? 



7-2 



A), the 



asymptotic form of self-similar solutions which approach 
the FRW solution for large z can be found by linearizing 
equations neglecting the V"^ term for accelerating expan- 
sion case {q > 1), while the l/V"^ term for deaccelerating 
expansion case (0 < g < 1). 



A. Accelerating expansion case 

For accelerating case (q > 1), the asymptotic behavior 
of linearized equations are given by 



C" = -il-q)C', 

2B" = -C" + -(k7 + 2(77 - 2 - k)A' 

7 



(4.4) 



+6{q-l)B' + {-q + K + KUj)C', (4.5) 

(4.6) 



^, ^ 7(K-l+tJAt) ^, 



(7-1)(2 + k) 



This system of linear differential equations has no solu- 
tion which vanish as z — > 00. It means that there are no 
non-trivial asymptotically flat FRW solutions in acceler- 
ating case. 



B. Deaccelerating expansion case 

For deaccelerating case {0 < q < 1) the asymptotic 
behavior of linearized equations arc given by 



-C" 



37- 



3 + 2cj 



(3r-|- 



3qK] C 



7 ■ 



-A 



+ i-2K+l-3q)C' - ""^^ ^h ', 



2B" = -C" + -(k7 + 297 - 2 - k)A' 

7 



-A' = 



7 



+6(g - 1)B' + {-q + K + KLu)C", 
-7 (3<z' - + 3qK^ (^C + 



-2(2 + k)B' - (k - 1 + ljk)C', 
which solutions is 



(4.7) 

(4.8) 
(4.9) 



C 



W + A^% (4.10) 

7 z** 7 z'^ 



s = l — 3q — K = —- 



2-7 


r , 2(37-5)1 
3(7-2)2 J 
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_'^+ 37(2-7)J 



(4.11) 



where Co,Ci and C2 are constants satisfying one con- 
straint equation 

- [uj{3-f^ - 87 + 8) + 4] Co + ^ [3a;7(7 - 2) + 37 - 8] Ci 



37-4 
37 



[3cj(7 - 2)2 - 2(37 - 5)] C2 = 0. 



The power s is positive when 

4 \ / 2(37-5) 



OJ - 



37(2-7) 



3(7-2)2 



< 0. 



(4.12) 



(4.13) 



When 7 = 4/3, s = —1/2 < 0. Since for deaccerelating 
expansion (0 < g < 1), we find 



to > 



(2-7)(37-2) 
1 



4 ^ 4 

> -, r tor 7 > -, 

37(2-7) 3' 



> 



■ 7 



4 ^ 4 



(4.14) 



the inequality (4.13) with (4.14) yields that the energy 
density of a perfect fluid is negative. Then the positivity 
of energy density requires s < 0, i.e., there is no non- 
trivial solution which approaches the flat FRW universe. 
We conclude that there are no non-trivial asymptotically 
flat FRW solutions in deaccelerating case. 
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V. DISCUSSION AND SUMMARY 

We have studied spherically symmetric self-similar so- 
lutions in Brans-Dicke theory which contain a perfect 
fluid with the equation of state p = (7 — l)fi (1 < 7 < 2). 
We have shown that when the energy density is positive, 
there are no non-trivial solutions which approach to the 
flat FRW spacetime at large z. It should be noted that 
although non-trivial solutions do exist, they are unphysi- 
cal solutions with negative energy density. It is similar to 
the result in general relativity. In general relativity, non- 
trivial solutions which approach to the flat FRW space- 
time at large z exist for 7 > 2, where the sound speed is 
faster than the speed of light ||] . 

The result in this paper suggests that the PBHs does 
not grow at the same rate as the size of the cosmological 
particle horizon in Brans-Dicke theory as in the case of 
general relativity. However, it could be possible to con- 
struct a non-trivial solution attached to the flat FRW 
solution at the point of = 1 or V^^ = 7 — 1 , where the 
field equations become singular. In general relativity, it 
was shown that the ingoing Vaidya null fluid solution 
can be attached to the flat FRW solution at the point of 
1/^ = 7 — 1 if the equation of state is p = /i though this 
is not a black hole solution j?] . 

If we allow deficit (surplus) angle in the present self- 
similar solutions, we can find a one-parameter family of 
solutions which approaches asymptotically to the "flat 
FRW spacetime", as Carr and Hawking found them in 
the case of general relativity [|| . 

In this paper, we have imposed the self-similarity on 
the whole spacetime, however it can be considered that 
the self-similarity is imposed only on the neighborhood 
of a black hole and the solution smoothly connects to 
the non-self-similar spacetime at large distance. Such a 
solution represents a black hole in the flat FRW universe 
which grows locally in the self-similar manner. To con- 
struct such a solution is an important future work. 
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